Transition of a mesoscopic bosonic gas into a Bose-Einstein condensate 
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The condensate number distribution during the transition of a dilute, weakly interacting gas 
oi N — 200 bosonic atoms into a Bose-Einstein condensate is modeled within number conserving 
master equation theory of Bose-Einstein condensation. Initial strong quantum fluctuations occuring 
during the exponential cycle of condensate growth reduce in a subsequent saturation stage, before 
the Bose gas finally relaxes towards the Gibbs-Boltzmann equilibrium. 
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I. INTRODUCTION 

Bose-Einstein condensates have turned into exquisite 
tools to study fundamental quantum phenomena on the 
micrometer scale, and a vast range of different physical 
scenarios have been realized experimentally with ultra- 
cold matter in the last decade, confirming the fundamen- 
tal importance and the broad applicational scope of Bose- 
Einstein condensation [l]. A microscopic, quantum dy- 
namical description of the gas' A/'-body state during the 
buildup of the condensed phase after a sudden switch of 
the gas temperature below the critical temperature ex- 
pected for Bose-Einstein condensation, however, is still 
one of the most striking theoretical topics of ultracold 
matter physics up to date. 

How can we model the quantum many-body dynam- 
ics during the transition of a dilute, weakly interacting 
gas of bosonic atoms into a Bose-Einstein condensate and 
link the resulting equilibrium statistics of the condensed 
phase to the roots of Bose statistics, i.e. to the statistics 
of quasi-ideal Bose gases, which was introduced almost 
one century ago by Bose and Einstein [2]? And under 
which conditions is the equilibrium statistics of a dilute 
and weakly interacting Bose-Einstein condensate of fixed 
particle number uniquely determined by the statistics of 
an ideal Gibbs-Boltzmann thermal gas - lacking any hys- 
teresis on the condensate formation process [3 ? 

Since an exact numerical solution to the full A^-body 
problem is out of the range of today's supercomputing 
facilities already for rather small atomic samples with 
a few hundreds of atoms, these questions reduce to find- 
ing appropriate and numerically accessible effective equa- 
tions for the quantum many-body dynamics of the Bose- 
Einstein phase transition. So far, extensive pioneering 
works [IHS] (and references cited therein) have led to ac- 
curate dynamic equations describing the evolution of the 
average condensate population for atomic gases with a 
few thousands to a few millions of atoms: the onset of 
Bose-Einstein condensation is marked by a spontaneously 
insetting exponential growth of the average ground state 
population, followed by a slow, subsequent saturation 
stage in which the condensate fraction converges towards 



an equilibrium value after suddenly cooling [7 the gas 
below its critical temperature. 

In this Article, the relaxation dynamics of the en- 
tire state of a mesoscopic Bose gas during the transition 
of exactly A" = 200 weakly interacting bosonic atoms 
into a Bose-Einstein condensate is reported for the first 
time within number-conserving master equation theory 
of Bose-Einstein condensation [131 US]- particular, 
focus is put on following the dynamics of the conden- 
sate number distribution during the Bose-Einstein con- 
densation process. Detailed theoretical knowledge on the 
statistics and dynamics of mesoscopic (A' ~ 200 — 1000 
atoms), weakly interacting quantum gases is not only 
of importance for state-of-the-art experiments [8^, but 
moreover contrasts fundamental postulates of classical, 
statistical mechanics to interacting, quantum degenerate 
many particle systems of finite size. The latter aspect 
arises from the neglect of number and energy fluctuations 
in standard thermodynamics, which is a reasonable as- 
sumption for classical, non-interacting gases in the ther- 
modynamic limit (with particle numbers as large as Avo- 
gadro's number, A^ ^ 10^^), whereas a vanishing impact 
of such quantum effects in the quantum degenerate limit 
is a priori no longer irrefragable. 

II. MASTER EQUATION THEORY 

In order to model Bose-Einstein condensation for small 
atomic gases, it is sufficient to numerically monitor the 
condensate and non-condensate number distribution dur- 
ing condensation, as will be shown in the following. For 
this purpose, we derive a master equation for said number 
distributions assuming a 3-dimensional harmonic trap- 
ping potential, with pn{N^^ t) = pn{N — A"^, t) the con- 
densate number distribution and A^^ the number of con- 
densed atoms, given a gas of A" particles therein. The 
particle number A" may be on the order of a few hun- 
dreds to a few millions of atoms, but is kept fixed during 
the final cycle of the condensation process. The deriva- 
tion of the master equation is lengthy and relies on sev- 
eral approximations [HI UBi , which can be summarized 
in the following way. We consider condensation onto one 
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FIG. 1. (color online) Schematics of microsopic many particle 
dynamics. The total number of atoms in the Bose gas is fixed 
to N. Thermal equilibration within the non-condensate is 
modeled by a heat reservoir which is at fixed temperature 
T. After suddendly cooling the gas, condensate and non- 
condensate undergo particle exchange due to the atomic two- 
body collisions until a steady state is reached. 



single-particle (system) mode 1^^), the condensate mode, 
whereas all orthogonal, non-condensate single-particle 
modes, {|^^),k = {kx^ky^kz)}^ are assumed to be lowly 
occupied (environment). Next, we employ the separa- 
tion of time scales in dilute atomic gases [Q [11] . The re- 
laxation dynamics in the high energetic, non-condensate 
part of the gas takes place on a time scale of two body col- 
lisions, Tcoiiision ^ ms, and is thus much faster than the 
time scale of condensate formation, tq ~ sec. Assum- 
ing this Markovian dynamics scenario which is basically 
justified by the rapid non-condensate thermalization, we 
formally map each Fock subspace of {N — N^) non- 
condensate particles onto a thermal mixture to model the 
rapid non-condensate rethermalization, and to care for 
particle number conservation. In particular, we replace 
interactions in the non-condensate thermal bath by for- 
mally coupling the non-condensate to an external heat 
bath of fixed temperature T. Finally, we use the Fock 
number representation, \N^^ {A/'/j}), with A"^, the number 
of condensate particles, and with A"^ = Xlfe/o ^fc' 
number of non-condensate particles. The thermodynam- 
ical scheme of our description is depicted in fig. [l] 

Since the particle number is conserved, the above in- 
gredients imply that the state of the Bose gas in Fock 
representation of the gas obeys merely classical correla- 
tions originating from particle number conservation. A 
Born-Markov ansatz [12] generalized for an A'-body state 
of a fixed particle number is thus given by: 



N 

pW(t)= pNiN^,t)\N^){N^\^p^{N-N^,T) . (1) 

N^=0 

Here, each Fock state |A^^)(A^^| is weighted with 
the condensate particle number distribution pN{N^^t)^ 



simultaneously capturing the distribution of non- 
condensate particles, as the particle number in the gas is 
conserved. Each state p±{N — N^^T) in eq. ([T]) denotes 
a thermal state projected onto the subspace of (A" — A^^) 
non-condensate particles, given that A"^ particles popu- 
late the condensate mode: 



'N-Nf> 



(2) 



with Z^{N-Ns, T) = exp[- ^£^0 /^^k^k]' the 

partition function of (A" — N^) indistinguishable parti- 
cles, in which the sum of the tuples {A"^^, A"^^, . . .} is 
taken such that the total atom number is conserved, 

Ek^k = {N- ^o)- The operator = Ek^o ^k4^k 
denotes the Hamiltonian of the non-condensate thermal 
vapor, /3 = (ksT)'^ the inverse temperature of the gas, 
and Qjv-ATg is a projector onto the Fock space of non- 
condensate number states with (A" — N^) particles. 

The fact that the N-body state in eq. ([2| is not a 
product state of a condensate and non-condensate den- 
sity matrix doesn't cause fundamental problems to derive 
a master equation for the reduced condensate density ma- 
trix. More explicit, as shown in Refs. [13l[T5], the state 
in eq. Q allows for the derivation of the master equa- 
tion for a three-dimensional harmonic trapping poten- 
tial without further approximations. Since the A'-body 
state is diagonal in Fock number represenation, we focus 
on the evolution of the diagonal elements characterized 
hyPNiN5,t) = {NolPoimo) = (iVolTr±/5(N)(t)|iVg) in 
eq. iQ. The corresponding master equation for the con- 
densate number distribution in a gas of exactly N atoms 
describes Bose-Einstein condensation as due to quantum 
jumps — j- ± 1 of the condensate particle number: 



dpN{NQ,t) 
dt 



[i+ {Nq, T) + e^(iVg, T)] Pm{Nq, t) 

+ ^+{NQ-l,T)pM{NQ-l,t) 
+ Q{NQ + l,T)pN{NQ + l,t) , 

(3) 

with a condensate feeding rate ^^(A'g,T) = 2{N^ + 
1)A+ (A"— A'g, T), and a condensate loss rate ^^(A/'g, T) = 
2N^XZ^{N — A'^,T). The single-particle transition rates 
A^(A' — A^^, T) are given in Eq. (|4|. 

Remarkably, even though eq. (|3| describes Bose- 
Einstein condensation realistically in terms of two-body 
atomic collisions ^ [H [13], it formally obeys the mas- 
ter equation for an ideal gas coupled to a thermal reser- 
voir J4]. In contrast to the non-interacting case, how- 
ever, the master eq. ([3| accounts for the real-time dy- 
namics of Bose-Einstein condensation in a dilute, weakly 
interacting gas. The master equation ([3| is valid up to 
order ag^^^ <C 1. Thus, even though we neglect terms 
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0{ag^^^) for the sake of simplicity, the time scale for 
condensate formation depends on both the type and the 
strength of the atomic two-body interactions. The con- 
densation time is quantified by the single-particle transi- 
tion rates A=^(A^ - N^,T) in eq. 



with /+(k, r, m, N±,T) given by 

/+ (k, r, m, N^,T) = {N^) {N^,T) {N^} {N^ , T) x 
x[{N^){N^,T) + l]\(,f\\ 



(4) 



(5) 



describing feeding processes of the condensate by one 
non-condensate particle, and with /_(k, 1, m, 7V^, T) 
given by 



/_(k,l,m,iVi,T) = [{Nk){N^,T) + 1] 



X [{Ni){N^,T) + 1] {N^){N^,T)\Cl{ 



mO|2 



(6) 



which quantifies losses of condensate particles towards 
the non-condensate single-particle modes originating 
from atomic interactions. 

Condensate feedings and losses with microscopic en- 
ergy balances + = ej^^ + eg ± fiA are quanti- 
fied in eqs. 6| proportional to the terms = 
/ dr ^j^(r)^|'(r)^^(r)^^(r), where = hkuj + ^(cj^ + 

ujy -\- ujz)/'^ with k = {kx^ky^kz) are single-particle ener- 
gies, and where A is the off-resonance of a given two- 
body process and F the energy uncertainty, in anal- 
ogy to the quantum optical case [12]. Thus, the S- 
distribution in eq. Q refiects and ensures energy con- 
servation on a certain width F arising from the finite de- 
cay time of non-condensate correlations [13j [15]. The 
{N^){N_\_,T) denote single-particle occupations, which 
are calculated with respect to each non-condensate ther- 
mal mixture of N± = {N — N^) particles in eq. ([2|, obey- 
ing Y.ki^k){N±,T) = {N - N^). The evaluation of F 
is beyond the scope of the present article, and wil be 
presented elsewhere [13]. 
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FIG. 2. (color online) Evolution of condensate particle num- 
ber distribution pN{N^,t) during the transition of a gas of 
N = 200 Rb^^ interacting atoms into a Bose-Einstein conden- 
sate of final temperature T/Tc = 0.4, obtained from eq. 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
time t [sec] 




time t [sec] 



FIG. 3. (color online) X-y projection of the distribution 
PN{N^,t) from fig. [2] (upper panel). Average condensate 
fraction ao{t) (lower left panel) and width ANo{t) in units 
{Nao{t)y^'^ of the distribution pN{N^,t) (lower right panel). 



III. DYNAMICS OF BOSE-EINSTEIN 
CONDENSATION 

The dynamics of the condensate number distribution 
during Bose-Einstein condensation is displayed in figs. [2] 
and [3] as a function of time t and condensate particle 
number A^^, after a sudden switch of the atomic cloud's 
temperature T below the ideal gas critical temperature 



Tc expected for Bose-Einstein condensation [T]. For nu- 
merical propagation of eq. (|3|, we consider a small Bose 
gas of = 200 weakly interacting ^^Rb atoms (with 
s-wave scattering length a = 5.4 nm [9] [ID]), prepared 
in a time- independent, three-dimensional harmonic trap- 
ping potential with trapping frequencies c«3 = {ujx = 
27r X 42.0 Hz, ujy = 27rx 42.0 Hz, cj^ = 27r x 120.0 Hz). De- 
spite the well-known S-shape behavior [4l, of the aver- 
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age condensate occupation number, the condensate num- 
ber distribution in particular gives rise of non-condensate 
number fluctuations during the condensation process [8] , 
since the total number of particles in the gas is fixed. 
Characteristic for a phase transition, we observe an ini- 
tial spread of the condensate number distribution (large 
number fiuctuations in the thermal vapor). While the 
condensate grows, these incipiently large fiuctuations re- 
duce until the reaching of a steady state. This steady 
state is in particular characterized in that the net parti- 
cle fiow between condensate and non-condensate is zero. 

Achieving Bose-Einstein condensation hence relies on 
the atomic interactions in order to bring the gas into 
equilibrium in the quantum degenerate regime, number 
fiuctuations are prominent during condensate formation 
as highlighted by figs. [2] and [3) and energy uncertainty 
(accounted for by the width T of the delta function in eq. 
Q) is present due to the finite decorrelation time be- 
tween condensate and non-condensate particles. Can it 
thus really be that - according to the fundamental laws 
of thermodynamics - a mesoscopic, weakly interacting 
Bose-Einstein condensate obeys the Gibbs-Boltzmann 
statistics? 



number of accessible microstates of the entire system, 
condensate and non-condensate. From eq. ([8|, it thus 
follows that the condensate statistics is independent on 
the specific strength of the two-body collisions once the 
steady state has been reached. This result is valid in 
the limit of dilute atomic gases and weak interactions, 
ag^^^ <C 1, as has been also verified numerically (innlud- 
ing the finite with F of condensate feeding and loss rates) . 
According to eqs. ([T]) and ([8|, the steady A/'-body state 
reached by the atomic collisions remarkably thus obeys 
the statistics of a thermal Boltzmann state of N bosonic, 
non-interacting particles [14]: 



N 



Z{N, T) 



Qn + 0{aQ^^^) , (9) 



with i-L = X^k^k^^^k' Hamiltonian of a non- 
interacting gas, Z{N^T) the partition function of N in- 
distinguishable particles, and Qat, the projector onto the 
Fock subspace of N particles. 



V. CONCLUSION 



IV. EQUILIBRIUM STATE 

To answer this question, let us recall that the steady 
state solution of the condensate number distribution in 
eq. ([3| entirely defines the equilibrium steady state of gas 
m eq. our framework, the key aspect that a Bose- 

Einstein condensate obeys the Boltzmann statistics lies 
in that quantum coherences in the gas rapidly decay on 
a certain time scale F~^ and in that the gas is sufficiently 
dilute ag^^^ <C 1. For sufficiently small energy uncertain- 
ties corresponding to the formal limit fSHT 0+, i.e. suf- 
ficiently fast non-condensate equilibration, phT <C 1, the 
rates A+(Ar- AT^, T) and A" {N-N^, T) obey the balance 
condition A+ [N - N^, T) = z{N - N^, T)\- {N - N^, T) , 
where 



^Z^{N-N^-1,T) 



Z^{N-N^,T) 



(7) 



In this case, the unique steady solution of eq. ([3| is de- 
termined by Pn{N^.T) ^ Uk'i' - hT)/^-{k,T). 
Using that A" (A" - A^^) X' {N - ^ 1), this steady 
state can be related to the balance condition in eq. ([7|), 
which leads to: 



^j,^Z^{N-Ns,T) 



Z{N, T) 



(8) 



where Z±{N — N^^T) is the partition function of the 
non-condensate gas in eq. ([2|, and where Z{N,T) is the 



In conclusion, we have monitored the condensate and 
non-condensate number distribution during the phase 
transition of a mesoscopic Bose gas into a Bose-Einstein 
condensate. We find that the steady state of the gas in 
the condensed phase is a Gibbs-Boltzmann thermal state, 
obeying the Bose-Einstein statistics of an ideal gas for 
sufficiently dilute and weakly interacting atomic gases. 
The reported scenarios were numerically reproducable for 
total particle numbers up to A^ = 10^ [15 . The impact 
of environmental decoherence sources onto the conden- 
sate formation process during condensate formation such 
as detailed study of deviations of the final steady state 
from the Gibbs-Boltzmann equilibrium will be of future 
interest and discussed elsewhere. 

I thank Andreas Buchleitner, Dominique Delande, 
Boris Fine, Cord Miiller, Alice Sinatra and Thomas 
Wellens for helpful discussions and comments during the 
development of the theory. Financial support from Hei- 
delberg University for the editing of this theory is 
acknowledged. 



[1] M. Anderson, J. Ensher, M. Matthews, C. Wiemann, 
and E. Cornell, Science 269, 198 (1995); K. Davis, M.- 
O. Mewes, M. Andrews, M.-O. Mewes, N. van Druten, 
D. Durfee, D. Kurn, and W. Ketterle, Phys. Rev. Lett. 
75, 22 (1995). 

[2] S. Bose, Z. Phys. 26, 178 (1928); A. Einstein, Sitzungs- 

ber. Kgl. Akad. Wiss. 1924, 261 (1924). 
[3] E. Altman, and E. Dernier, Nature 449, 296 (2007). 



5 



[4] C. W. Gardiner et. al, Phys. Rev. A 55, 4 (1997); 

D. Jaksch et. al, Phys. Rev. A 56, 575 (1997); C. W. 

Gardiner et. al, Phys. Rev. Lett. 79, 1793 (1997); C. W. 

Gardiner et. al Phys. Rev. A 58, 536 (1998); D. Jaksch 

et. al, Phys. Rev. A 58, 1450 (1998); D. Jaksch et. al, 

Phys. Rev. A 61, 033601 (2000); C. W. Gardiner et. al, 

Phys. Rev. Lett. 81, 5266 (1998); M. D. Lee et. al, Phys. 

Rev. A 62, 033606 (2000); M. J. Davis et. al, Phys. Rev. 

A 62, 063608 (2000). 
[5] H. T. C. Stoof, Phys. Rev. Lett. 66, 3148 (1991). 
[6] R. Walser, J. Williams, J. Cooper, and M. Holland Phys. 

Rev. A 59, 3878 (1999). 
[7] S. Ritter, E. Ottl, T. Donner, T. Bourdel, M. Kohl, and 

T. Esslinger Phys. Rev. Lett. 98, 090402 (2007). 
[8] R. Gati, B. Hemmerling, J. Foiling, M. Albiez, 

M. K. Oberthaler, Phys. Rev. Lett. 96, 130404 (2006). 
[9] Ch. Buggle, J. Leonard, W. von Kitzhng, and J. T. Wal- 

raven, Phys. Rev. Lett. 93, 173202 (2004). 
[10] S. Stringari and L. Pitaevskii, Bose-Einstein condensa- 



tion (Oxford Science Publications, 2003). 

[11] H. -J. Miesner, D. M. Stamper-Kurn, M. R. Andrews, 
D. S. Durfee, S. Inouye, and W. Ketterle, Science 273, 
1005-1007 (1998). 

[12] C. Cohen-Tannoudji, J. Dupont-Roc, and G. Grynberg, 
Processus d 'interaction entre photons et atomes (Savoirs 
Actuels, Editions du CNRS, 1988). 

[13] A. Schelle, T. Wellens, D. Delande, and A. Buchleitner, 
Phys. Rev. A 83, 013615 (2011) 

[14] VI. V. Kocharovsky, Vi. V. Kocharovsky, M. Holthaus, 
C. H. Raymond Ooi, A. A. Svidzinsky, W. Ketterle, and 
M. O. Scully, Advances in Atomic, Molecular and Optical 
Physics 291, v.53 (2006). 

[15] A. Schelle Formation of a dilute Bose-Einstein con- 
densate: Number- conserving master equation theory of 
Bose-Einstein condensation (Siidwestdeutscher Verlag 
fiir Hochschulschriften, ISBN 978-3-8381-1289-3, 2009). 



